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Abstract. By using a linearized non- vacuum late time solution in Brans-Dicke 
cosmology we account for the seventy five percent dark energy contribution but not 
O |' for approximately twenty-three percent dark matter contribution to the present day 

■ energy density of the universe. 
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iy-j ■ Our universe seems, according to the present-day evidence, to be spatially flat and 

to possess a non vanishing cosmological constant [U I2J- For a flat matter dominated 
universe, cosmological measurements [2j imply that the fraction of the contribution 
of the cosmological constant A to present energy density of the universe is ^a ~ 0.75. 
In standard cosmology ^a would be induced by a cosmological constant which is a 
dimensionful parameter with units of (length) -2 . From the point of view of classical 
general relativity, there is no preferred choice for what the length scale defined by A 
might be. Particle physics, however, gives a different point of view to the issue. The 
cosmological constant turns out to be a measure of the energy density of the vacuum 
^ ■ and although we can not calculate the vacuum energy with any confidence, this allows 
us to consider the scales of various contributions to the cosmological constant. The 
energy scale of the constituent (s) of A which in Planck units is approximated to 10~ 123 
is problematic since it is lower than the normal energy scale predicted by most particle 
physics models. To solve this problem, a dynamical A [4] in the form of scalar field 
with some self interacting potential [E] can be considered and its slowly varying energy 
density induces a cosmological constant. This idea called " quintessence" [3] is similar 
to the inflationary phase of the early universe with the difference that it evolves at a 
much lower energy density scale. The energy density of this field has to evolve in such 
a way that it becomes comparable with the mass density fraction now. This type 
of specific evolution, better known as " cosmic coincidence " m problem, needs several 
constraints and fine tuning of parameters for the potential used to model quintessence 
with minimally coupled scalar field. To solve the cosmic coincidence problem, a new 
form of quintessence field called the "tracker field" [TT] has been proposed. Such kind 
of quintessence field is mainly based on an equation of motion with a solution for such 
that for a wide range of initial conditions the equation of motion converge to the same 
solution. This type of solution is also called an 'attractor like' solution. There are a 
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number of quintessence models proposed. Most of these involve minimally coupled scalar 
field with different potentials dominating over the kinetic energy of the field. Purely 
exponential jSj-[Z| and inverse power law [Ej-[H] potentials have been extensively studied 
for quintessence fields to solve the cosmic coincidence problem. However the fact that 
the energy density is not enough to make up for the missing part of the cosmological 
constant or that the p/p = 7 value found for the equation of state of quintessence is 
not in good agreement with the observed results makes such an explanation unlikely. 
The investigation of alternative models in which the equation of state parameter 7 of 
the cosmological constant evolves with time has been proposed due to the conceptual 
difficulties associated with a cosmological constant [T2]-[TK]. 

There have been quite a few attempts for treating this problem with non-minimal 
coupled scalar fields. Studies made by Bartolo et al [01], Bertolami et al Ritis et al 
[TH] have found tracking solutions in scalar tensor theories with different types of power 
law potential. In another work, Sen et al [19 have found the potential relevant to power 
law expansion in Brans-Dicke (BD) cosmology and Ank et al |2Hj have shown that (BD) 
theory of gravity with the standard mass term potential \m 2 (j) 2 is a natural model to 
explain the rapid primordial inflation and the observed slow late-time inflation. 

In this paper we show that a linearized non-vacuum solution about the stable 
cosmological vacuum solution with flat space-like section is capable of explaining how 
the Hubble parameter evolves with the scale size of the universe a(t). In this framework, 
we also show that the standard Friedmann equation changes into a form in which the 
power of the scale size term with Qm is corrected by an amount 1/u 

Hy n A + n M (^) 3+TS (i) 



Ho ) u 

where u is the Brans-Dicke parameter with uj 3> 1 |2*T) \'2'2\. Subsequently, under such a 
linearized solution, we point out that only a very small part of the dark matter can be 
accommodated into the contribution of the Brans-Dicke scalar field. 

In the context of (BD) theory [231 with self interacting potential and matter field, 
the action in the canonical form is given by 



S = J d 4 x yfg 



^(t> 2 R+\ <T drf dA - \m 2 <p 2 + L M 
z>ui 2 2 



(2) 



In particular we may expect that <fi is spatially uniform, but varies slowly with time. 
The nonminimal coupling term <p 2 R where R is the Ricci scalar, replaces with the 
Einstein-Hilbert term jjr-R in such a way that G~^ = ^4> 2 where G e ff is the effective 
gravitational constant as long as the dynamical scalar field <fi varies slowly. In units 
where c = H = 1, we define Planck-length, L p , in such a way that L 2 p (pQ = u/2n where 
(f)Q is the present value of the scalar field (p. Thus, the dimension of the scalar field 
is chosen to be Lp 1 so that G e ff has a dimension Lp. The signs of the non- minimal 

coupling term and the kinetic energy term are properly adopted to (H ) metric 

signature. The Lagrangian of the scalar field, in addition to non-minimal coupling term 
and the kinetic term, is composed of a potential which consists of only a standard mass 
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term. Lm, on the other hand, is the matter part of the Lagrangian which in accordance 
with the weak equivalence principle is decoupled from as has been assumed in the 
original (BD) theory. Excluding 0, as the matter field, we consider a classical perfect 
fluid with the energy-momentum tensor Tjf = diag (p, —p, —p, —p) where p is the energy 
density, p is the pressure. 

The gravitational field equations derived from the variation of the action (j2J) with 
respect to Robertson- Walker metric is 

3 a (a 2 , k \ 1 12 1 2 ,2 3 a ■ 
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where is the curvature parameter with k = —1, 0, 1 corresponding to open, flat, 
closed universes respectively and a (t) is the scale factor of the universe (dot denotes 
4). Since in the standard theory of gravitation, the total energy density p is assumed 
to be composed of p = p^ + pm where p^ is the energy density of the universe due to the 
cosmological constant which in modern terminology is called as "dark energy", the right 
hand sides of (jHlllJ) are adopted to the matter energy density term p M instead of p and 
Pm instead of p where M denotes everything except the field. The main reason behind 
doing such an organization is that whether if the terms on the left-hand side of (jHJ) can 
accommodate a contribution to due to what is called dark matter. In addition, the right 
hand side of the equation (jSJ) is set to be zero according to the assumption imposed 
on the matter Lagrangian Lm being independent of the scalar field 0. By defining the 
fractional rate of change of as F (a) = 0/0 and the Hubble parameter as H (a) = a/ a, 
we rewrite the left hand-side of the field equations (jHHSl) in terms of H(a), F(a) and 
their derivatives with respect to the scale size of an universe a (prime denotes 4-) 



,2 2uj 2 ^ jj 1 2u 2 / 4oA p M 



H*+l— + l] F 2 +-HF+— [HH + HF\ + — -—m 2 = — ^-(7) 
V3 3y 3 3V /3a 2 3 \ 3 / ^ U 

H 2 -^F^uHF + a (^~^HF] +±-^ m 2 = . ( 8 ) 
3 I 2 3 J 2a 2 3 K ' 

Solving the field equations (jUHHI) for the closed universe stable-vacuum solution, 
we get = cf) p e Fpt and a (t) = a* w 1/y/um where Fp ~ 0.7m is the primordial 
fractional rate of change of 0, a* is the constant size of this static universe, P is the 
constant value of the field as t —>■ and m is the mass of the scalar field 0. The point 
to note here is that the closed universe (k = 1) vacuum solution becomes important 
for the primordial universe since homogeneity of the universe only makes sense if a 
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closed universe undergoes big-bang. Under these closed universe vacuum solutions, we 
have shown [20] how the presence of radiation changes the behavior of the universe 
compared to these stable solutions. Putting the scalar field vacuum solution ~ e Fpt 
in (J3j) and imposing initial conditions as the big-bang time and the corresponding time 
dependence of the scale size of the universe, we have found the solution for the scale 
size of the universe in the form of 



where c is an integration constant and H p is the primordial Hubble parameter with 
H p = ujF p 0.7ujm for u ^> 1. This solution provides importance for the following 
reasons: 

(1) It is a natural solution. Namely, it does not need any "special" equation of state 
for the matter. It is just deduced from the theory by putting the stable-empty universe 
solution (ft ~ e Fpt into the equation (jSJ) 20J . 

(2) If one examines this inflationary solution concerning as t — > and as t > 0, it 
is seen that (jUJ) is both consistent with a (t) ~ \ft as t — ► and also with primordial 
rapid inflation described by a (t) ~ e Hpt for uj 3> 1. 

(3) We have also checked that for u ^> 1 and in the limit as t — > 0, if one 
substitutes 4> ~ e Fpt and a ~ \/t into (JHH^) then the equation of state p = l/3p is 
satisfied automatically as expected in the radiation dominated epoch of the standard 
Einstein cosmology. 

In the light of this encouraging result obtained by using the instability caused by 
the nonvacuum in the closed stable vacuum solution in explaining the rapid primordial 
inflation, we will show that a linearized non-vacuum solution about the flat stable 
vacuum solution can also be powerful in explaining the slow late time expansion. 
Since the universe becomes (approximately) flat in late times, we ignore the curvature 
parameter k/a 2 as a (t) increases with the expansion of the universe. Under these 
considerations, in analogy with the assumption we use in explaining rapid primordial 
inflation, we first propose a = e Hoct and cf) = e Foat and put into (jUKl) and search for 
a zeroth order stable vacuum (empty except the field) solution. Hoo, are the 
constants to be determined named as the late time Hubble parameter and the fractional 
rate of change of in the late time regime respectively. We have the following coupled 
equations for and F^; 




(9) 




^luF^ + 2H 00 F c 



oo 




(10) 






oo 




(12) 



and get the solution; 




(13) 
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Foe = — x m « — = m (14) 

where the approximations are again for u ^> 1. Thus, in our one hand, we have had an 
exact zeroth order stable- vacuum solution as, 

a = e Hact (15) 
,e F ^ (16) 



where i?oo ~ O.Sy/ujm, Fao ^m,, 0oo is a constant. 

Then, after finding such a zeroth order exact stable solution, the question that 
stimulates us, similar to the primordial regime analysis, is that how the presence 
of matter affects this flat stable vacuum solution (JT51 [Tfij) . To understand such a 
perturbation phenomenon, we impose the following linearized first order non-vacuum 
solution for H = a/a and F = </>/</> which includes first order perturbation functions of 
h(a) and f(a) in addition to the constant terms and which appear in the flat 
stable vacuum solution (|T5| ITH|l respectively. 

H = H 00 + h(a) (17) 

F = F 00 + f(a). (18) 

Since solving the field equations (jBHSJ) exactly for a(t) and under the condition p = 
is hard enough, we put our imposed solution (fTTl HH1 ) into the modified field equations 
(J3EJ) for p = and neglect higher terms in h(a), f(a) then we get h(a) and f(a) for all 
uj in the form of, 

Ha) = Cl g„ (!*)**-( ^- s ) <" + V'r ± | (i9) 



o/ yH^alJ (uj + 2)(2uj + 3) U 

3^+4 



/(a) = C7 2 ^o(^)' :++1 + 



. , (20) 
(cj + 2)(2o; + 3) V a 7 v 7 

where ao is present size of the universe and H is the present Hubble parameter. C\ and 
C*2 are, on the other hand, dimensionless integration constants. Since letting uj — > oo 
has a special meaning in the sense that the Brans-Dicke scalar tensor theory matches 
with standard Einstein theory under such limit, we display the linearized solution (JT71 
ITS]) in the following form as uj — > oo, 



a ^ 3+1 



a 



F = F oo + C 2 H (^) . (22) 



Hence, putting the solution (|21jl in the standard Friedmann equation, 



-i =iu+n«^ +«M^j (23) 
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which is used for fitting Hubble parameter to the measured density parameters of 
universe in such a way that Qa+ ^_r+ — 1 and using the present observational 
results on density parameters ~ 0.75, Qm ~ 0.25, Qr ~ 3J, we get G\ ~ 0.15 
and ifoo = 0.86iJ so that « H 00 /2uj rs (0.43/a;) # which provides |C 2 | < 0.43/a;. 
Namely, the first term in a linearized solution (|22jl is much greater than the second term. 
The curvature density parameter Qr, on the other hand, is found to be in accordance 
with the recent measurements since the term (l/if ao) 2 ~ &r — [3]. 

To compare (jHJ), with standard Friedmann-Lamaitre cosmology, we put the 
linearized solution (|2*T1 l2*2*j) into this equation and transfer all terms except for H 2 = 
(a/ a) 2 to the right hand side. Neglecting the 1/a 2 term 
4a; 

H 2 = —(p A + p M + p D ). (24) 

Noting that, in the late time regime, ~ a 1//2w is approximately constant as a changes 
we identify the terms which do not explicitly depend on a with p\ and terms which 
depend on a as a -3 with the dark matter energy density po so that 

Pd = (C^H^ 2 - ^H^ 2 - ^H.F^C^f (25) 

2lo Zoo a 

PA = \Flo<f ~ ^#°°Foo0 2 + \ m2( t> 2 ' ( 26 ) 

Using the recent observational results on density parameters of the universe (Qa 
= Pa/ Po — 0.75, = Pd/ Po — 0.23) where p is the present measured energy density 
of the universe and the relations F^ m 11^/200 ~ (0A3/u)H and « 0.8y/um as 
to — > oo, we fit (jSniEEJ) to the ratio ^a/^d — 75/23 and determine the \C<z\ integration 
constant to be | C 2 1 ~ 0.20 which is inconsistent with the requirement |C 2 1 <C 0.43/a> 
imposed by the theory. 

In conclusion, the first remarkable feature of this work that a linearized non- vacuum 
solution (|2*T|) about the stable cosmological vacuum solution (JT3J) with flat (k = 0) space- 
like section is capable of explaining how the Hubble parameter H = a/a evolves with 
the scale size of the universe a(t). 

The second remarkable feature of this theory is that by fitting the linearized 
solutions (|2*T1 l2*2*jl of the theory to the recent observations j^j, the late-time Hubble 
parameter = 0.86i?o and the fractional rate of change of <fr i n the late time regime 
Foo = (0.43/a;) H are successfully predicted in terms of today's observational measured 
vale of Hubble parameter H . 

Another important prediction we note from this theory is that for a fixed Hq, since 
Foo ~ (0.43/a;) H , F^ may not attain a large value because of its inverse dependence on 
uj which is measured to be, according to the recent observational data, as u > 10 4 3> 1 
[2*T| This is the reason why F^ can not let the scalar field <fi = $ OQ e F ° ct to blow 
up rapidly so that pa (ESI), the energy density of the universe due to the cosmological 
constant, can grow slowly and reasonably. Hence we strictly agree on that this theory is 
successful in explaining the dark energy though the Brans-Dicke scalar field <fi can not 
account for dark matter. 
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The last remarkable feature of this theory is that it enables us to estimate some 
dimensionful parameters displayed in the theory. Using the relation H^ w 0.86ifo ~ 
0.8y/um and the restriction on u, we may estimate m for a fixed H as 

m < IQ~ 2 H (27) 

where the present value of Hubble constant H = 720 ± 8 km s _1 Mpc _1 [24 . Using 
appropriate conversion relation in relativistic units, present Hubble parameter is found 
to be H « lO-^m- 1 « 2 x 1(T 42 GeV. 

By using the relation G~h = (2ir /u)(j) 2 , the time variation in Newtonian gravitation 

constant can be written as | G/G \= 2F = 2 \J>/4>) so that we can estimate the time 
variation in Newtonian constant G, Fp, F^ for the primordial and the late time epochs 



considered in this theory as 

|^|p = 2Fp^(1.5/v^)# (28) 

17^00 = 2^00^ H (29) 

F P < 7 x 1(T 3 H (30) 

< 43 x 10- 6 H . (31) 

The Hubble parameters, on the other hand, are given by 

Hp w 0.7muj > 70 H (32) 

#oo w 0.8^^ w 0.86iJo- (33) 
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